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SUMMARY 


It is shown that the correlation of fluctuating static pressure (in 
an incompressible and homogeneous turbulence) with any fluctuating quan- 
tity in the flow field, can be expressed in terms of the correlation of 
the same quantity with two or more components of the velocity. 

The correlations of pressure with itself and of pressure with two 
velocity components are investigated in detail for the case of isotropic 
turbulence. These correlations can be expressed in terms of correlations 
involving two velocity components at a point and two velocity components 
at another point. A postulated relation between the fourth-order and 
second-order correlations is investigated. This relation is satisfied, 
for example, if the joint probability density of the four components of 
velocity is Gaussian. The consequences of this relation are compared 
with the measurements of the fourth-order correlations. 

The root-mean-square pressure and pressure gradients are computed 
from second-order correlation for a range of turbulence Reynolds num- 
bers. Since the pressure gradient is related to diffusion of marked 
particles from a source, the computed pressure-gradient level is com- 
pared with that calculated from a set of diffusion measurements. 

The triple correlation equation and plausible hypotheses relating 
higher order correlations with second-order correlation are examined for 
the possibility of getting a determinant set of equations for isotropic 
turbulence. 


INTRODUCTION 


Although in isotropic turbulence the dynamical (correlation) equa- 
tion does not contain any pressure correlations, the correlation of pres- 
sure with itself is of interest in turbulent diffusion and in the study 
of sound generation by turbulence (for very low Mach numbers of turbulence 
velocities when there is little effect of sound on the turbulent field, 
that is, only a slight amount of energy is drained away from the turbulent 
energy in the form of sound waves). The static pressure fluctuations are 
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also Int ima tely connected with the onset of cavitation in turbulent flow 
of liquids. 

The pressure-gradient fluctuation, in connection with diffusion in 
isotropic turbulence, was first considered by Taylor (ref. l). He postu- 


lated that — 



?ldx 


is of the order u n 


o/^Un 




independent of Reynolds 


number of turbulence"'". Heisenberg (ref. 2) later derived, from detailed 
spectral considerations, an expression for as follows: 

( 1 / p 2)(^/ Slx l) 2 “ Ul ( Sui / 5x l) _ where = u 2 /\(dui/dx : x ) 2 is the 

(0.13) 2 R a / n ' / 

Reynolds number of turbulence and the symbols are defined in the symbol 
list. 


Obukhov (ref. 3 ) derived an expression for the correlation of pres- 
sure with itself for distances between two points whose spacing is within 
Kolmogoroff ’ s "inertial subrange." Batchelor (ref. 4) and limber (ref. 5 ) 
have considered the pressure-pressure and pres sure -and-two-velocity- 
component correlations in isotropic turbulence. Chandrasekhar (ref. 6) 
has considered these correlations for the isotropic turbulence in 
magnetohydrodynamics. All these correlations can be expressed in terms 
of the correlations involving two conponents of velocity at one point 
and two conponents at another point. These investigators have from the 
very beginning assumed a simple relation between the fourth-order and 
the second-order correlations. In order to compare the primary result 
with experiment it is necessary to have the results in terms of the 
fourth-order correlation.^ 

Next to the simplest case of istropic turbulence is homogeneous 
axisymmetric turbulence, in which all statistical properties are symmetric 
about a particular axis, instead of being spherically symmetric as in the 
case of isotropic turbulence. This introduces one new element: The 

energies in different conponents of the velocity are not the same and 
there is a transfer of energy from one component of velocity to another. 


In the present paper an overbar indicates mean value or statistical 
average and an underbar indicates a vector quantity. 

^The main results were derived before the publication of refer- 
ences 4 to 6, but publication was delayed until the conpletion of the 
experiments reported here. 
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This transfer is essentially due to the correlation "between velocity 
and pressure or pressure gradient (ref. 7)* This is in addition to the 
transfer of energy from big eddies to small eddies ; both of these 
transfers are expressible in terms of the nonlinear terms of the equa- 
tions of motion and these can be expressed in terms of the triple -order 
correlations. Axisymmetric turbulence is the only case of homogeneous 
turbulence that has been theoretically studied and even this has been 
largely limited to general tensor forms for the correlations and equa- 
tions governing these correlations. The physics and solution of the 
problem remain virtually untouched. In the general case of homogeneous 
turbulence the pressure-velocity correlation is undoubtedly of interest. 

Before some special correlations are taken up, correlations involving 
static pressure in general will be considered and it will be shown that 
these correlations can be expressed in terms of correlations involving 
more than two velocity components. This is convenient from the theor- 
etical point of view, since it may then be possible to express higher- 
order correlations in terms of second-order correlations on the basis of 
some plausible hypotheses, such as those used by other investigators. It 
may be also noted that Heisenberg’s expression for the spectral transfer 
term (ref. 2) is a hypothesis for the triple correlation in terms of the 
second-order correlations. The fact that correlations involving static 
pressure can be expressed in terms of correlations involving more than 
two components of velocity is also convenient from an experimental point 
of view, since at present there exists no technique for the measurement 
of static pressure fluctuations. On the other hand, the standard hot- 
wire technique can be extended to correlations involving more than 
two velocity components. 

This work was sponsored and financially supported by the National 
Advisory Committee for Aeronautics. The author should like to thank 
Drs. S. Corrsin, F. H. Clauser, and L. S. G. Kovasznay for their helpful 
suggestions, and to thank the following people for their help in measure- 
ments and data processing: Mr. A. L. Kistler, Miss Patricia Clarken, 

Miss Patricia O'Brien, and Miss Ingeborg Busemann. 


SYMBOIS 


Al,A2 

a,b,c,d 

a i> b i 

C 


dimensional constants 
arbitrary unit vectors 

direction cosines of two arbitrary directions at x 
constant 
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absolute constant (eq. (28)) 
velocity of light 

direction cosines of two arbitrary directions at x 

electric-field strength 
output of hot-wire 

output of hot-wires set at x and x', respectively 

P( r ).»g( r ) correlations used by Von Karman and Howarth 
H magnetic -field strength 

h = (p./ 4 np)H 

j current in conducting fluid in motion 

K universal constant 

L = v5A e V4 

M mesh spacing 

n integer 

P = (l/p)(p - P) + ( 1 / 2 ) (|h| 2 - jh| 2 ) 

P' 

P 

q' 

R 

r 1 , 2,3 
R A . 
r 


value of P at x' 

static pressure fluctuations 

any quantity at another point x' 

scalar correlation 
correlations 

Reynolds number of turbulence 

distance between two points x and x 1 > (^i^i) 


c i' d i 


I 

e 

e l> e 2 
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3 1. . .n 


t 

U 

u 


u 


v (s) 

v 

V 1 

w 

x i 


set of random variables 
time 

velocity in x-direction 
velocity component 

components of instantaneous fluctuating velocity, 
or u, v, w 

volume 

velocity in direction perpendicular to hot-wire 
root-mean-square velocity in y-direction 

coordinates of point x 

£ 


u l> 


U2, 


y 

A = 
5 


variable of integration 
Sx^Sx-l 

Kronecker delta 


A 

\ 

v 

£ 


energy dissipation per unit mass of fluid, (l/2) (d^u^/dt 
Euler ian microscale, 

Lagrangian microscale 

permeability 
kinematic viscosity 
component of displacement vector 
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p density 

p charge 

0 conductivity 

r time difference 

^(S^,S 2 ,S^,Si ( .) joint probability density of S-j_ . . . S^ 

characteristic function of <f>", tp, t2> t^, and % 

are independent variables 

(_) vector 

( ) statistical average or mean 

Subscripts: 
e excess 

Z,m,n 

PjOj^s 

P 


free indices characterizing general vector component; each can 
take value of 1, 2, or 3 corresponding to component in x-, 
j- } and z-directions, respectively 

directions (see fig. l) 

free indices 

refers to pressure 
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GENERAL RELATIONS 


Consider a statistically homogeneous and incompressible turbulence 
with no mean velocity (i.e., a "box turbulence"). The equations of 
motion for an arbitrary point x are for momentum 

K/at) ♦ u 3 SJ - ^ + v (l) 

and for continuity 


cki^/dx-^ = 0 


(la) 


where the symbols have the usual meaning and repeated index means sum- 
mation over the index. Taking the divergence of equation (l) and making 
lose of equation (la), the following well-known relation results: 

1 a 2 P _ (2) 

^ dx^ dx^ cbc^ 8x^ 


The fact that the pressure in an incompressible viscous fluid satisfies 
Poisson's equations means that it is not a primary variable, since the 
pressure can be expressed in terms of the velocity derivatives by using 

the well-known solution for Poisson's equation. 5 

Multiplying equation (2) by q', any quantity at another point x* , 
and taking a statistical average. 


1 , ^p 

— q« £ — 

P cbc^ Bxj_ 


a 2 - 


-<3/ 


UiUj 


dx^ Sxj 


The solution of the equation A n )6 = r(x) where A = 

is discussed in reference 8 . In the following discussion the solutions 
for n = 1 (Poisson's equation) and n = 2, in three dimensions, are 
needed. Here 0 and r are arbitrary functions of 2 and n is an 
integer. 
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Ma kin g use of the fact that the turbulence is statistically homogeneous 
and interchanging the order of differentiation and averaging. 


1 _ ^UjU^' 


where p(x)q' (x* ) = pq’(|) and x 1 = x + _£. An overbar denotes statis- 
tical average. Equation ( 3 ) is just the Poisson equation for pq' and 
its solution (ref. 8) can he written immediately as 


1 

P 


pq'(t) 


1 _ 

r ° 

8 UiU^q' (jr) dV(y) 

J 

Syi |y - jj 


(4) 


where dV(y) is a volume element. Equation ( 3 ) in one form or another 
is the basis of the recent investigations reported in references 3 to 6 
which have been mentioned earlier. Next particular cases of equation (3) 
are considered. 


Case (l): q’ = u k ’ 


When q* = u k ’. 


2 

1 r = _l_ 5 Uj^jV dV (il) 

P PU]£ 4rtJ 8y ± y-J[| 


The quantity pu k ' is- zero in isotropic turbulence because it is a first- 

order isotropic tensor, but it is different .from zero in the general case 
of homogeneous turbulence. Its role in the axisymmetric turbulence is 

discussed in reference 7- lu axisymmetric turbulence pu k ' gives the 

transfer of energy from one component of velocity to another. As in 
isotropic turbulence there is transfer of energy from the large eddies to 

small eddies which is given directly by UjUjU k '. This transfer and the 
transfer of energy from one component to another can both be expressed 
in terms of u^u .u^. 1 . 



X 


MCA TIT 3116 


9 


Case (2): 


q' = - 


d 2 p' 




q uk'^z' 

^k'^z' 


Substituting q’ = — 


d 2 p' 


P ck.£ 1 8x.j_ 1 


5 w 

^k’^z' 


into equation ( 3 ), inter- 


changing the order of differentiation, and averaging. 


1^. 

d pp> 


= 5 

p 2 a s± a E± at, aq aq as a lk a Ej 


(6) 


Making use of the fact that the elementary solution of the bi-Iaplacian 
in three dimensions is r/2 where r = is the distance between 

the point of observation and the point of integration (ref. 8), 


ijjTij— i- r±wvv_ | y .s| av(y) 

d 2 8n J c^. ay. a>v — ~ 


(7) 


Case ( 3 ): q' = u k' u z' 

When q' = u^'u^ ', 


1 

P 


PVV 


1_ r^VVVV mte) 


( 8 ) 


Case (4): Correlations Involving Pressure in an 

Incompressible, Highly Conducting Fluid 

/ 

The correlations involving pressure in an incompressible, highly 
conducting fluid have been investigated in reference 6. The current Jt 
in a conducting fluid in motion is given by 


= tf(cE + (iu x H) + (Pei i/ c ) 



10 


MCA TN 3116 


where 0 is the conductivity, E is the electric -fie Id strength, u 
is velocity vector, H is the magnetic-field strength, p is the perme- 
ability, c is the velocity of light, and p e is the excess charge. 

The quantity P e H/ c represents the convection current. As 0 — > 00 it 
may he assumed that 


E « -pu x H/c 


otherwise the current will become large . The energy in the electric 

field is of the order |n| 2 /c^ of " t * e ener Sy in the magnetic field, 
therefore it can be negelected. In this approximation only the inter- 
action between the two fields u and H needs to be considered, and 
the equations of motion are (ref. 6) 


^i | d 

8t dx,, 


( u i u j - h i h j) ■-? + ?N 2 ) + v Au i 


(9) 


1 + T-( h ± h 5 - UiUj) = A Ah ± 


dt 2bc r 


( 10 ) 


where h = \h and A = -A itpc. Taking the divergence of equation (9) 
- \hxpj- 4 

and making use of the condition of incompressibility (la). 


2 

cgp _ 3 - ^3) 


cbc^ dx^ 


dxi dxj 


( 11 ) 


where 


P = i(p _ p) + i(|h [ 2 - |h| 2 ) 


Chandrasekhar (ref. 6) has considered three correlations Pu-^'u^', Ph-^'lij 1 
nnfl pp 1 . it can be seen from equations (7) and (8) that these correlations 
can be expressed in terms of u^UjU^'u 


l> h i h jW' and u i u J h k ,h z'- 
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The quantities Puj^Uj' and pu^Uj' have s imi lar tensor forms and they 
satisfy similar differential equations (equations (ll) and (2)). Also 

PP' and pp' satisfy similar differential equations, so that it is not 
necessary to calculate separately all these expressions. When the 

expressions for pp' and pu-^'^' have been calculated the corre- 
sponding expressions for other correlations can be written at once. 
Furthermore, pp’ and pu^'^' are not independent. If q 1 = p 1 in 
equation (3 ) 

1 jjji jvg 

In all these cases (except case (l)) the correlations involving 
pressure can be expressed in terms of correlations between two velocity 
components (or two magnetic -field components ) at one point and s imilar 
quantities at another point. Some investigators have assumed that the 
joint probability density of these four quantities is Gaussian (refs. 4 
and 9 )- If Sp, S£, S^, and are any such variables then it can 

be shown that 


SiS^Si, — SpSj + SpSj^ + SpSg (12) 

and in particular 

“AW “ v7 + VV AV + \’“ I ' <«> 

h i h jW - h lV + h lV h j h k' + h l h ,1 VV (W) 

u i u A'Y (15) 

(The opportunity has been taken to correct the expression corresponding to 
equation ( 15 ) in reference 6 where the last term u_jUj h-^'h^’ = u^h-j^S^S^ 
(no summation) is missingV 
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The hypothesis of Gaussian joint probability density makes all odd 
moments (e. g. , triple correlations) zero. However, the above relations 
are only integral conditions on the joint probability density and the 
hypothesis of Gaussian probability density is strictly not necessary 
for their validity. Without making any such hypothesis about the prob- 
ability density, equation (12) is taken as a plausible hypothesis. 

The joint probability density )6(Sp, S 2 ,S 3 ,S]j.) is the Fourier trans- 
form of the characteristic function t 2 , t^, t^) 

*(^ 2 ,^,%) = JJJJH S^Sg^S^e 131 ^ dS x dS 2 dS 5 dS 4 

"|( s k s J t k t l) 

= e ‘ L 


where (f> is Gaussian and \|r is a function of independent variables tp 
to %. The moments of all orders are determined from the behavior of 

the characteristic function near the origin. As has been pointed out by 
Chandrasekhar (ref. 6), if 


*(tl,t 2 ,t 3 ,t 4 ) 


^SkS^kh)^ + Cp^tptqtr) 


is taken as the characteristic function, then there is the same relation 
between second- and fourth -order moment as for a Gaussian probability 
density, but the third-order moment will not, in general, be zero. In 
view of this it can be asserted that the hypothesis of equation (12) 
does not impose any restriction on the triple correlations. 

Cases (2) and ( 3 ) will be investigated in detail for the simple 
case of isotropic turbulence and the results compared with experimental 
observations . 


FOURTH-ORDER VELOCITY CORRELATIONS IN ISOTROPIC TURBULENCE 


The fourth-order correlation 



enters prominently in the 


expressions for the correlations involving static pressure fluctuations. 
It is a fourth-order isotropic tensor and its form can be derived using 
some simple results of the invariant theory (see appendix). Thus: 
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u i u jV u 2 


22, \ nn, . 22, \ , n2 

R,,(r) + R (r) - 2R (r) - Ij-R , 
22 nrr ' nn v ' n2 


MjM2 + ^m( r ^ij^k2 + 


t[o-> - + s ii 5 ^) + 


22 

nn 


(r) - 


c 2*ni l (r) - «r) - 

2r _ 


+ W ±3 )*-^ 

+ ? 3 5 k 8 ll + + "J 


Since u^u-ju^'u^ ' is not solenoidal, the continuity equation gives no 

relation between the five scalars defining u^UjU^'u^ 1 . However, it 

gives relations between derivatives of u^u^u^'u^' and other fourth- 

order correlations involving derivatives of velocity. Some representative 
correlations are .shown in figure 1. 

Under the assumption of equation {12.), 


*22^ 


«*0 


R n J(r) 

n2 


= 2 


+ R 


22 


M 

2 l§ (r 3 2+B ^ : 




> 


^ 1 


R™ = R R 
nn nn mm 


(IT) 




These relations and the expression for u^u^u^'u^' for the case k = 2 = 
were first given by Millionshtchikov (ref. 9 ). The second-order corre- 
lations R n (r) and R^(r) (essentially the Von Karman-Howarth g(r) 
n L- 
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f(r)) are connected by the continuity equation of incompressible 
flow (ref. 10). 

£<*> “ *!<*> + I | R l ( r) (18) 

With the corresponding inverse relation, 

R*(r) = \ f l£(y)y dy (l8a) 

0 


The fourth-order correlation u^UjU^u^ * involving three velocity 

components at one point and one at another is not directly related to 
the pressure correlations that have been considered. It enters explicitly 
in the equation for the propagation of triple -order correlation which will 

be considered later. For isotropic turbulence, the correlation u^uju-^u^ ' 

has the following form (see the appendix): 

• u i U ,jW = + Cn - 3E Un ~ 5 R ?nn)MA^ + 

f E mm( S lJ S kz + + 8 iZ 8 kj) + 

^Zn - 3 Cn)(*i S Al + WjZ + 5 d S Az) + 

(*?nn - 5 C.)(WiJ + Wik * Mz 5 Jk) (19) 

The various correlations sire shown in figure 1. Here u-jUjU^u^' is a 
solenoidal tensor and the four correlations characterizing it are not 
independent. Taking the divergence of u^UjU^u^’ with respect to the 
last index and equating it to zero. 



> 


( 20 ) 
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Because of these two relations there are only two independent correlations 
characterizing the tensor u^uju^u^'. It is observed in passing that 

UiUiU k u z ' is a solenoidal second-order tensor. This gives a relation 

between the four functions characterizing u^u^u^n^ ' , which, however, is 

not independent of those written above and serves as a check on the alge- 
braic errors. If the hypothesis of equation (12) is satisfied. 


E m * 


Cn = 


r“ = R l7 R n 
ZZn 11 n 


R? = R R* 
Znn nn Z 


-A 




( 21 ) 


STATIC-PRESSURE CORRELATION IN ISOTROPIC TURBULENCE 


For the isotropic case 


) = j] + b d 3 ( ) _ 1 d- 4 [r( 

^b 

and the bi-Iaplacian equation (6) simplifies to 


r dr 5 


r ,4 
dr 


1 a rpp'(r) _ d UjUjVu^ 


p 2 r 


Ctr 1 * Slj & s k 


( 22 ) 


After four repeated integrations, for any fixed value of t, corresponding 
to equation ("j), 


PP 


Jt 

- 1 r m * o U.U.Uv'U, 1 

* = — / y(y - r ) 5 1 J k 1 

6r J r Sy . Sy , Sy, c 


\ 


ay 


(23) 


It was assumed that pp' > 0 fast enough as r — > 00. Substituting 

equation (l6) in equation (23) and carrying out some algebraic simplifi- 
cations and integrations by parts. 
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PP'(r) = R u (r) " 




(24) 


The pressure gradient is of special interest since it is related to 
the diffusion of marked fluid particles from a fixed source (see the 
next section, "Pressure Gradient and Diffusion Measurements"). Now, 


where 


^_\ 2 = lim fr (* i ) M * 1 ') 

dx-J r — >0 dxi ^1' 


x^' = x^ + r 


Interchanging the order of integration and averaging. 


/a? \ 2 _ 

^PP’(r) 

V ^1/ 

_ dr 2 


r=0 


Substituting for pp’ from equation (24), 




* A\ - < 


nl \ dy 

y5 

(25) 


Tn order for the above integrals to converge it is necessary that the 

series expansion of (r™ 1 + R*J - 2R U - 4R^? ) at the origin begin with 

\ nn ll nn n IJ 

terms of order r^j therefore (r™ 1 + R^ - 2R^ - 4R^) = 0. Use 

<3^2 \ nn ll nn 

has been made of this fact in deriving expression (25). Under the assump- 
tion of equation (12), equations (24) and (25) become 
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and 


-4 PP'(r) 


■/>#) 




— , 2 


dy 


dy 




2 


_d_ 

dy 


B^(y) 


dy 

y 


(24a) 


(25a) 


For vanishing turbulence Reynolds number. 



where 



and 



The correlation R^(r) can be calculated neglecting the inertia terms 

L 

of the equations of motion (ref. 11), 

R*(r) = u-l 2 exp (-r 2 /2A 2 ) 

Substituting this result in equations (24a) and (25a), 


pp'(r) = P 2 ^ 2 ] exp (-r 2 A 2 ) = R* (r) 

and 



(26) 


(27) 


Op the basis of Kolmogoroff 1 s hypothesis of local isotropy, for large 
Reynolds numbers and values of r with the "inertial subrange" (ref. 12) 
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(u - u 1 ) 2 = 2 


U* - Ej(r)] = Cl (ve) l/2 (E) 2/5 


(28) 


where u is the velocity along the displacement vector, = x' - x, 


P 1 du i u i 

r = 1^1^, e = — — is the energy dissipation per unit mass of the 

fluid, the length L = v^^e'* 7 ^, and C-]_ is an absolute constant. The 
quantities e and 7\ are connected by Taylor's decay equation (ref. l): 


1 = -WA 2 

2 dt 1 ‘ 


Hence for values of r within the inertial subrange 


M.- 2 

2p P 


P^ - PP' 


)-/ 

n 


r r% 7 - 

s'(y) 
dy l 


- r “n, , n 2 

yay - r U 


^ . (29) 

y 


The integrals in equation (29) can be evaluated approximately by using 
equation (28). After some transformation the result is 


I p 2 (p 2 - pp' ) - \ «(|) 


,v 3 


or 


i(p-P') = 


(u - U' ) C 


(50) 


This relation was first given by Obukhov (ref. 3) and equations (24a ) 
and (25a) were given by Batchelor (ref. 4). 


PRESSURE GRADIENT AND DIFFUSION MEASUREMENTS 

In reference (13) the "Iagrangian" microscale was related to 

the acceleration in isotropic turbulence without the neglect of viscosity. 
A brief resume follows: 



NACA IN 3116 


19 



(31) 


where v' is the root-mean-square velocity in the y-direction. Sq uar ing 
and averaging the v-component equation of motion. 


(to? f _ 1 /^P\ 2 

Vdt / p2\Sy/ 


+ v 


■(Av) + 


— .^E 
P by 


Av 


(32) 


The term 


(!-)= 


0 because of isotropy. All terms in the expansion 


of equation (32) can he expressed in terms of the derivatives of R^(r) 
by making use of the relations given by Von Karman and Hcwarth (ref. 10 ). 


Thus (Av)^ = — 

3 


a^(r) 


dr* 


Townsend (ref. 14) has made measurements 


r=0 


of 


Sr* 1 


and his result is 


r=0 




R z ( r ) 


— 1 r=0 


_ 35 (v 1 f ( 30 . 

3 *4 W 


0. 2R-j 


(33) 


For the isotropic turbulence 
formulated by Taylor in ref. 


cl(v' ) 2 

dt 


= -lOv 


l). Thus 



(the "decay equation" 



( 34 ) 
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Substituting equations (32), 


(33), and 


(3^) in equation ( 31 ) 


k 






H. _ 

_ 2 R-x 

r A A I 


-1 


( 35 ) 


Heisenberg (ref. 2) has calculated an expression for pressure gradient 
using the hypothesis that various Fourier components of the velocity 
field, are statistically independent. As far as the fourth-order corre- 
lations are concerned the consequences of this hypothesis can be shown 
to he equivalent to equation (12) (see ref. 4). His result is 



3K 

78.5 


R-> 


( 36 ) 


Instead of Heisenberg's original extrapolation formula for the stationary 
velocity spectrum, Chandrasekhar's solution (ref. l6) was used to calcu- 
late equation ( 36 ). Here K is a "universal" constant which has to be 
determined from experimental data. Heisenberg found from turbulence 
decay that K = 0. 85 . lee (ref. 17 ) gives its value as 0. 13 based on 


^This is the first step toward the establishment of a possible rela- 
tion between Iagrangian and Eulerian correlations. If nondecaying iso- 
tropic turbulence is considered then the Lagrangian correlation 


u i(t) u i(t + t) = B^t) = u x 2 + 


(-D 



T 2 * 1 

(2n)J 


Each term of this series can be expressed in terms of the Eulerian corre- 
lations by using equations of motion. The first nonconstant term involves 
pressure or fourth-order velocity correlations ► The next term will involve 
still higher order correlations. It may be possible to relate higher order 
correlations with second-order correlations by making a further hypothesis 
about the joint probability density of the velocities at two points. The 
first step of this computation has just been carried out. The next step 
will involve many more calculations. This is probably not the best method 
of attack of the problem of relating the Lagrangian and the Eulerian corre- 
lations. Frehkiel has related Lagrangian and Eulerian microscales under 
some assumptions (ref. 15 ). 
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fba i\ 

skewness of the probahility density of — as measured by Townsend 

A^i 1 

(ref. 14). Froudman (ref. 18) has calculated double and triple corre- 
lations from Chandrasekhar's solution of Heisenberg's self -preserving 
spectrum. He has estimated its value as 0.^5 from comparison of theo- 
retical correlations with measured correlation. Heisenberg's results 
are only approximate and the value of K depends on the range of spec- 
trum (or correlation) which is made to fit the experimental data. The 
value of K = 0A5 has been taken as a compromise value. 


PEESSURE -VELOC ITY CORRELATIONS pu k 'u x 


Equation (8) relates the correlation pu^'u^ ' with the fourth-order 
correlations 


1 ^pVV _ d^ujW 

p ^i 


(37) 


Here pu^'u^' 
same form as. 


is a second-order isotropic tensor; 


for example, the tensor u^u^ ' 


therefore it has the 



Z Z yj-p 

The correlations Rp (r ) and Rp (r ) 

tuting the expression for u^uju^'u^ ' 
of equation (37)> 


are shown in figure 1. Substi- 
(eq. (l6)) in the right-hand side 
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5u 1 U A' U Z 

aii^j 


3r 2 


Is]] - S U ) * i |-(e“ + 4s“ - SR 

\ ZZ nn ) r Sr \ I'- 


ll 

'nn 


nZ mm\ p / nn ZZ __ZZ -,^-rP-l 
12E nZ + R nn) + ^pl^Rnn + R ^Z “ ^ Rnn ” 12511 1 ~ Rrm / 


£k^Z 


A 1_ r f2R n J + 5^ - R m U 1 *(seP + 1*; 

2r ^2 \ nj nn nn/ r orl m 


«- «)+£$£♦< - 


5 kZ 


= -^(r) - t 2 (r)5 kI 


(39) 


where and. tyg are functions of r which are defined by the above 

equation. Substituting equations ( 38 ) and (39) in equation (37) > 


1 _d_ 
r dr 


1 

I? 


|^(i'] r3 ) Mi + Ipr + ^sWi * -♦1 - Mti 


Equating coefficients of the above equation. 


_d_ 

dr 


3 3?( p i r5 ) = -|V r) 


After integrating twice. 


p,r- 


- -i /> - r?K 


dy 

y 


m 


end 


APo = « 


1 u rp2 2Pl 

“ -*2 - -r 


2 r Sr 2 
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After integrating twice. 


•• ■/. (*> * 5 


3 


[7 - — ) dy 


However, it is more convenient to calculate + 3Pgj the correlation 

between static and dynamic pressure. It follows directly from equa- 
tions (37) and (39) that 


a(Pi + 3p 2 ) - i ^(r, ( Pl + 3p 2 )j = -(*! + 31-2) 


After integrating twice. 


p l + 3p 2 


'/ (*1 + 3 *2)( y - t ) ** 


(41) 


Substituting for i|r-^ in equation (40) and integrating by parts. 


Pi(r) = 


ll nn 


5r 5J 0 
12r 5 


fC + + K! 


5FT 1 - 

nn 


O" * - fc * B}J - < + <}) ^ (42) 


and from equation Oil), 


P 1 + 5p 2 


_ 1 


% P u k' u k 


- - R U - + 3 (“?f - 2 f (> 


K 1 ” - h” - h" + b”) 32 (43) 

nn nn ll nn/ y 
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From equation ( 38 ) and equations (42) and (43 ), 

nn Z 2 ( o\ 2 1 r 1 "/ nn 2 Z __ nZ CTi m 

Rp ~ ~^nn " \ U 1 ) x / (^^nn + ^®2Z + ^®nZ " + 

* VyP J o 

dy - ^jT“(C * - a» - O | - 

f/ r "(C-^-^-C)f <»> 

and 

+ (“?f + ^5 /'(< + - 

5 R ™ 1 _ 4R U )y 2 <3y + | f (r 1 ^ + R l ] - 2R^ - 
nn nn J 5 J \ nn ll nn 

^ r 


Using the hypothesis of equation (12), equations (43), (44), and (45) 
become 


1 

P 


pW 


j.i 

Pp + 


nn P° 

-± 


i •:« 2 


y <*y 


^ - 





(43a) 

(44a) 

(45a) 
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Comparing equations (24a) and (43a) it is seen that 

\ P 2 = "%> u k u k (^6) 

This result is a consequence of the equations of motion and of equa- 
tion (12). 

As in the case of pressure-pressure correlation, an explicit rela- 
tion for pressure-velocity correlations for the limiting case of small 
and large Reynolds numbers can be derived (see eq. (2 6 ) and the related 
text ) . For R^ > 0, 

R^(r) = u x 2 exp (-r 2 /2A 2 ) 

Substituting this in equation (43a), 



For R^ — > and, on the basis of Kolmogoroff ' s hypothesis of local 

similarity, for values of r within the inertial subrange (see eqs. ( 28 ), 
( 29 ), and ( 30 ) and the related text) 

(u - u' ) 2 = 2 J^Rj (r ) + u 2 J = c(ve) 1 / 2 ^ 2 ^ 5 (48) 

Substituting this result in equation (43a), 

|(pu k u k - pu k 'u k ' ) = JL C 2 vg(£) 2 / 3 (49) 

or using equation (48) 

^(p u k u k - pVV ) = 12 



( 50 ) 
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or using equation (30 ) 

i(pu k u k - pu k 'u k ' ) = ^(p2 ~ PP' ) (51) 


Equations (bjis.), (k4a), and (k5a) and (k9), (50), and ( 51 ) were first 
given by Limber (ref. 5)* 


TRIPLE-CORRELATION EQUATION IN ISOTROPIC TURBULENCE 
Consider the equation of motion at the point x and time t 


duj duj^j = _ 1 ap_ 

at Sxj p Stx^ 


Multiplying this equation by 



and taking the statistical average. 


Uv'U 


k * aT 


dui aUiUjUj,'^ 


cbc4 


1 Hv 

P SbCjL 


+ v Au^'u^ 


(52) 


In order to express 



in terms of the triple correlation 


u^u it is necessary to proceed as follows: Consider the equations 

of motion at x' and t: 


dy + sy v = _i v_ + 

at ax s 1 p ax k 1 


v A'u k ' 


+ dyv i ap' 

at ax s ’ p ax^ 


t 


+ v A'Uj ' 
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Multiply the first equation by 1 and the second equation by u k ' 
and add the two equations: 

Vv 4. iHAilV - . iL. <*' + v ^1, + 


dt 


SXs' 


= r ax. 




k UA Z 

v(v A'V + V A'uj ') 

Multiplying the last equation by u-j_ and taking the statistical average. 


du,, 


k'V ^k U Z ' u s ,u i 


i,-, «„ 


U i St 


dXc 


= -? U ± r i ' l ^ + U *'^) + 


v(u z ' A'u k ' + U k ' A , u z ')u i 


(53) 


to 


It is not possible to express u ^ 1 A'u^.'^ in terms of u^'u^’u^ or 
dp' 

express Uj ^ ^ in terms of Uj'p'i^ without introducing new 
^k' 

correlations. The attempt to express u k 'u^ ' in terms of u k 'u i 'u.^ 

has resulted in introducing two new correlations, each of which is sole- 
noidal third-order isotropic tensor characterized by a single scalar. 
Adding equations (52) and (53) gives the equation for the propagation 

Of UjUfc.'Uj' 


^iW + ^ u i u ,i u k tu z' + ^ u i u k' u Z ,u s' _ _ 1 ^P u k' u Z ' _ 
dt dxj dx s ' P dxi 


p(v + V t^rVi + v A Vk'V + V (Y + V a’Uj*)^ 


(5^) 
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Millionshtchikov has derived an equation relating third- and fourth- 
order correlations (ref. 9 ). His result is 


duk'V^i doiUjW 


8t 




1 dpyy 


v Au^y 



Millionshtchikov used this equation together with the hypothesis of 
equation (12) to derive an expression for the second-order correlation 
which is valid when the inertia terms are small. As a first approxi- 
mation he calculated the double correlation from the Karman-Howarth 
equation by neglecting the triple correlation, which is valid as 
R-^ 0. As a second approximation he used his equation for the triple 

correlation and the hypothesis of equation (12) to express the triple 
correlation in terms of double correlation. In his calculation he neg- 
lected the correlation pu^'u^’} however, this omission can be rectified 

since pu^'u^' can be expressed in terms of the second-order correlation 

by using equations of motion and the hypothesis of equation (12). The 
error in his equation for the triple correlation is more serious. When 
this is corrected a rather complicated equation for the triple corre- 
lation results (eq. (53) above). 

The purpose in deriving an equation for the triple correlation was 
to express the triple correlation in terms of correlations which can be 
expressed in terms of double correlation by use of a plausible hypothesis. 
If this were possible there would be (with the Karman-Hcwarth equation) 
two equations for two unknowns, that is, the two scalar functions char- 
acterizing double and triple correlation tensors. This has been only 

partly successful, however, since u^u-ju^y , u^u^u^'u^ ' , and pu^'u^' 

can only be expressed in terms of double correlations by using the 
hypothesis of equation (12) and the equation of motion. The two remaining 
correlations appearing in the equation for the third-order correlation 
cannot be expressed in terms of known correlations without making further 
assumptions. Some of these difficulties can be overcome if correlations 
can be considered that involve quantities at three different points and 
later make two of the points coincide. 

Previously It was believed that higher-order correlation equations 
could not be investigated because of pressure-velocity correlations 
entering in these equations . Most of the pressure-velocity correlations 
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can be expressed in terms of correlations involving velocity components 
at two points, and others, in terms of velocity correlations at more 

than two points. Consider, for example, the correlation pu^Uj 7 which 


is related to u H Ui'. The latter quantity enters in the triple- 

dx k 1 J 


correlation equation. Now, pu.u.' = lim pu. "u, ' and it is seen from 

3 x"-^x 3 

equation (3 ) that 


d p(x)u^(x" )uj (x 1 ) S u^.CxOu^x^Cx'^UjCx' ) 

= -P 


cSbti cbc^ 


^k ^2 


so that by using three point correlations, it is possible to express all 
the pressure-velocity correlations entering in the triple -correlation 
equation in terms of velocity correlation. It may be possible to relate 


u^u^u^'Uj" with double correlations u k a i ', u^'u^", u k u j"' 80 
forth, using, for instance, the hypothesis of equation (12). 


It is of interest to investigate what limitations and errors, if 
any, are introduced in the dynamics of turbulence by the use of the 
hypothesis of equation (12). If the fourth-order correlations could be 
simply related to the triple correlations it would be possible to inves- 
tigate theoretically some of the consequences of the hypothesis of equa- 
tion (12). 


MEASUREMENT OF SECOND- AND FOURTH- ORDER CORRELATIONS 
Equipment and Technique 


The measurements of the second- and the fourth-order correlations 
were made in a 2- by 2-foot wind tunnel at 48 mesh lengths downstream 
from a 1-inch-square mesh grid made from l/4-inch circular rods. The 
electronic equipment used is described in reference 19 . The new equip- 
ment used here is an electronic squaring circuit which has been developed 
by Kovasznay (ref. 20). The circuits consist of pairs of rectifiers with 
series resistors acting as a full -wave rectifier. All pairs are in par- 
allel; however, each pair is biased more than the preceding, so that as 
the input voltage increases more and more of these pairs of rectifiers 
conduct. The bias voltage and the series resistors are chosen so that 
the total rectified current of all the conducting pairs is proportional 
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to the square of the instantaneous voltage. The squaring circuit responds 
ins tantaneously within the limitations of the capacity effects in the 
diodes. Calibration with a harmonic signal showed that its response is 
good up to 70 kilocycles. 

The schematic diagram of the equipment used in measuring the corre- 
lations is shown in figure 2. The preamplified hot-wire signal is fed 
to a push-pull power amplifier which drives the squaring circuit. The 
output of the squaring circuit is passed through a microammeter and a 
thermocouple. Because of the inertia of the moving coil of the ammeter 
it is sensitive only to the average current or mean square of the input 
voltage. The output of thermocouple, due to lag, is proportional to the 
average square of the current and hence it is proportional to the average 
fourth -power of the input voltage. This arrangement measures simulta- 
neously the average square and average fourth power of the input voltage 
and this enables one to measure simultaneously the fourth-order and the 
corresponding second-order correlations. 

The squaring circuit has hardly any measurable error. However, the 
thermocouple, because of radiation and other losses, is not a perfect 
squaring device, and the error depends on the probability density of the 
signal. The complete fourth-power circuit (the squaring circuit and the 
thermocouple) was calibrated with a "noise" generator which gives a 
signal having Gaussian probability density (see fig. 3)- In view of the 
fact that in most of the turbulence measurements approximately Gaussian 
signals are dealt with, the above calibration was used to correct all the 
measured data. The error in the corrected results is expected to be 
within ±5 percent. If ej is the output of the hot-wire set at x 

and e2 is the output of the hot-wire set at x' and furthermore if 
the two wires are set perpendicular to the mean flow, then 

e l = a l u 

e 2 = a 2 u ' 

where u is the fluctuating velocity at x in the direction of the mean 
flow and u 1 is the corresponding quantity at x'j then 


( e l + e 2) 2 ~ ( e l ~ e 2) 2 _ e l e 2 = uu' 

k\ e-L 2 e ^ 2 ) ej 2 e 2 2 V u 2 (u 1 ) 2 
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( e l + 62 )^ + ( e l " e 2 ) J+ ' “ 2& 2^ ei 2 e 2 2 u 2 (u' ) 2 

and 


( e l * e 2p ~ ( e l ~ e 2)^~ _ e l^ e 2 + e l e 2^ u(u 1 P 

8(1? ^) lA (5T + ^f) 1/2 J?^ A (^4 + ^) 1/2 ? 


if u(u' = u^u 


so that in order to measure the fourth-order corre- 


lations u 2 (u ' ) 2 and u^u and the second -order correlation uu' it 
is only necessary to measure the average square and average fourth power 
of (a) the outputs of two hot-wires, (t>) their sums, and (c) their 
difference. 


The correlations R**(r) and R* (: r ) were measured by moving 
one hot-wire with respect to the other in the direction of the mean flow. 

The correlations R ^ (r ) and ®^ nn ( r ) were measured by moving 

one hot-wire with respect to the other in a direction perpendicular to 
the mean flow. 


The combination ( R + R + 2R + 4R ) was measured by moving 
\ 44 nn l L n 4 / 

one hot-wire with respect to the other along a line mating an angle of 
45° with the mean motion (see fig. 8 ). 

The combination (r?, + R™ 1 + 2R™j 1 - ^R n »') was measured by using 
\ 4 4 mi 4 4 n 4 / 

two wires perpendicular to each other and each inc lin ed to the mean 

flow. When one wire is on top of the other the combination of wires is 
essentially an x-type hot-wire which is sensitive of v-conponent of the 
fluctuating velocity. 


The correlation R was measured by using one x-type hot-wire 

sensitive to v and a single wire sensitive to u and moving apart the 
two probes in the direction of v. 
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The correlation was measured by using an x-type hot-wire sensi- 

tive to w and a single wire sensitive to u and moving the two probes 
apart in the direction of v. 

The length of the hot-wire was quite small (the ratio of the micro- 
scale ~k to the length being about 10 ), so that no appreciable length 
correction was necessary. 


Results 


The measurement of various correlations was made in isotropic tur- 
bulence at R^ = 60 . The correlations R^ (r ) and R^(r ) were meas- 
ured (see fig. 4) and they satisfy reasonably well the relation 

-R“(r) 


which arises from the continuity equation of the incompressible flow and 
the condition of isotropy (ref. 10 ). This served as a check on the iso- 
tropy of the -turbulence. 


The correlations 


ll , . nn. . mm. . 

(r ), R nTi (r ), and R nn ( r ) and two independent 


combinations 


'll nn ll , n l 

R,, + R + 2R + 4R , 
^ ll nn nn n l 


^ and 


ll nn ll n V 

R,, + R + 2R - hR , 
ll nn nn nZ, 


were measured. The measurements are compared with the hypothesis of 
equation (12) in figures 5 to 9* For large displacement of the points 

the correlation between u^ and disappears and u-j^Uj ' ^ 

tends to a constant value u _^ 2 ^ u j ' ^ = (n^) Tor isotropic turbulence. 
If the probability density of u^ is assumed to be Gaussian then 
n2 


V 


- ?( 


* 1 2 ) 


and therefore for large displacement of the points 



a dashed line in figures (5), ( 6 ), and ( 8 ). While the fourth-order 
correlation was measured the corresponding second-order correlation was 

also measured simultaneously. The correlation R™? (fig. 10 ) was also 

L L 

measured; this is not independent of the above five correlations. S in ce 
the hypothesis of equation ( 12 ) gives it a constant value it is of interest 
to see if the measured value deviates from this constant. Within the 
experimental accuracy, all the fourth-order correlations are close to 
those conputed from corresponding second-order correlations using equa- 
tion ( 12 ). 
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X zi 

The correlations R^^(r) and R mm (r) characterizing the quad- 
ruple correlation u^uju^uj 1 were measured. These are compared with 

those computed from corresponding second-order correlations using equa- 
tion (12) in figures 11 and 12. The hypothesis of equation (12 ) is 

satisfied for R^,, and R n within experimental scatter. 

Ill nnn 

CALCULATION OF ROOT-MEAN-SQUARE PRESSURE AND PRESSURE GRADIENT 
FROM VELOCITY CORRELATIONS AND DIFFUSION MEASUREMENTS 


Let R-i , Ro. and R* denote the correlations R^ and R 1111 and 

0 ll nn 

the combination j~/R^ + R™ + 2R^ + 4R^^.> res P ec, tiv e ly. This set of 

correlations, among other possible sets, suffices to determine the pres- 
sure correlation. In terms of the above notation equation (24) becomes 


4> PP'(r) = R x (r) - (u-] 2 ) 2 + 2^ (3R2 + R x - 4R 5 ) - 


^(R-L + Rs - 2R 5 ) 


r L_ 

£y 
y 


where R2 was measured for both positive and negative values of r 

(see fig. 6) and the point r = 0 was determined from the fact that it 
is an even function of r. The correlations R-[_ and R^ were measured 

for positive values of r and there is some uncertainty (0.025 inch) in 
the determination of the point r = 0 for these two correlations. This 
uncertainty, the experimental scatter, and the fact that small differ- 
ences between relatively large quantities are required to compute the 
pressure correlation have made the results very uncertain. 

Two widely different sets of curves were drawn for each measured 
correlation, taking into account the experimental scatter and the uncer- 
tainty in the determination of the point r = 0. Two pressure corre- 
lations were computed from these two different sets of fourth-order 
correlations. The result of these computations and those from the 
second-order correlation are compared in table I and figure (13). 
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It appears that for the computation of pressure correlation one 
should measure dir ectly the differences ( 3^2 + Rp ~ 4R^) aad 
(R-L + R 2 - 2Rj) instead of the individual correlations. This involves 

the difficult problem of measuring small correlation between relatively 
large quantities. In view of the fact that the hypothesis of equa- 
tion ( 12 ) is approximately satisfied, the double correlation can be used 
to compute root-mean-square pressure and pressure gradient. However, 
one unfortunate circumstance must not be overlooked. It is possible that 
the hypothesis of equation ( 12 ) is approximately satisfied but the root- 
mean-square pressure and the pressure gradient computed from using this 
hypothesis are still in error because the root-mean-square pressure and 
pressure gradient depend on the differences of various fourth-order 
correlations . 


Equations (24a) and ( 25 a) express the desired quantities in terms 

of the longitudinal second-order correlation . The correlations R^ 

i 0 

and R n are connected by equations (l 8 ) and (l 8 a). It is necessary to 


n 


express 


— R^ in terms of R 11 . Since 
dr l n 


dr 


14s 1 




therefore 


_d_ 

dr 


= ~ f ^ 

0 r ? '-'O ay 


Using this result, equations (24a) and ( 25 a) become 
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Root-mean-square pressure and pressure gradient were computed from meas- 
ured values of R^(r) for various Reynolds numbers. The results of the 

computation are shown in figures 1^ and 15 . The data for the two highest 
Reynolds numbers were taken from reference 21. The data in reference 21 
are given in terms of grid Reynolds numbers. Turbulence Reynolds num- 
bers R^ and A/M (where M is the mesh spacing) were est ima ted from 

grid Reynolds numbers by assuming that U^/u-^ t - -'- in the initial 

period where U is the mean velocity. 

Equation (35) relates the Lagrangian microscale with the root- 

mean-square pressure gradients. Measurements of have been made by 

Simmons (reported in ref. l) and Collis (ref. 22 .), and an extensive set 
of measurements is given in reference (13). The experimental data 
(table II ) are used to compute the pressure gradient (eq. (35)) and. the 
results are compared with Heisenberg's analysis (eq. ( 36 )) in figure 15* 


CONCLUDING REMARKS 


It is of theoretical interest and experimentally convenient that 
correlations involving static pressure fluctuations can be expressed in 
terms of higher velocity correlations j the latter in turn can be related 
to second-order velocity correlations by using some plausible hypothesis. 
The fourth-order correlation, which enters prominently in the corre- 
lations pp' and pu^'u^ ' , may be related to the second-order corre- 
lation by using the hypothesis of equation (12). Experiments lend sup- 
port to this hypothesis. Since the differences of quadruple correlations 

are involved in the expression for pp ' , a slight deviation of the quad- 
ruple correlations from those computed using double correlations can lead 

to considerable error in pp'. Further experimental improvement involves 
the difficult problem of the measurement of small correlation between 
relatively large quantities. 

Ci 

The experimental determination of root-mean-square pressure gradient 
from diffusion measurements is also not very accurate because of the 
inherent double differentiation of the experimental curves involved in 
this technique. The hypothesis of equation (12) relates the Lagrangian 
microscale to Eulerian microscale. It may be possible to relate the 
Lagrangian and Eulerian correlations for larger values of the independent 
variables by use of the equations of motion and higher-order correlations. 

The dynamic equation relating quadruple correlations to triple corre- 
lations is worth further investigation, firstly, to see if one can possibly 
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get a closed system of equations for the dynamic of isotropic turbulence 
hy use of the hypothesis of equation (l 2), and, secondly, to see the 
limitations and errors introduced hy the use of the above hypothesis. In 
this respect it is noted that experiments lend support to the hypothesis 

of equation (12) for UjUjU^'Uj ', the quadruple correlation involving 


two components at one point and two at another, and also for UjUjUjjUj ' , 

the quadruple correlation involving three components at one point and 
one at another. If this were strictly true then the quantity 


( u i - V)V ( U 1 - u i')5 


would he constant and equal to a numerical 


value of 3 (the value for Gaussian joint distribution ) independent of 
position, where Uj_ and Uj_' are velocity components at two separate 


points and perpendicular to the displacement vector. The quantity 



has been measured (fig. l6) and its value 


deviates from 3 tor small displacement of the points and in this region 


it approximately equals 



/cJujA 2 


\^*2/ 


The deviation of 


(u - u' ) y [(u - u' )^j from the numerical value of 3 for small dis- 
placement of points shows that the hypothesis of equation (12) is not 
satisfied for small eddies; however, the maximum error is only about 
20 percent. Since the differences of correlations are involved in the 

quantity (u - u' ) VEu-u.fJ small deviations in the values of 
UfUjU^’u^ ' and u^uju^u^' from that computed from the hypothesis of 

equation (12) show up prominently in the values for (u - u' )^/jT u - u')^j 
while ^these small deviations are hardly noticeable in the values for 

u i u jW “ d u i u jVV- 


Batchelor (ref. ^ ) has discussed the correlation u^u^u^'u^ 1 in 
connection with pressure fluctuations. He has presented Stewart's meas- 


urements of 
density of 


(u - u' ) 

Uf, u 3 . 



U£ ' and u-j 1 


to show that the joint probability 
is Gaussian. As far as the pressure 
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fluctuations are concerned it is only necessary to assume that u^uju^'u^' 

satisfies the hypothesis of equation (12) which is less restrictive than 
the hypothesis of Gaussian joint probability density of u^, u ^ , u^'. 


and u^ ' . 


Even if 


(u - u* )^j (u - u' ) 2 


is equal to a numerical value 


of 3 this does not prove that u^ Uj, u k ', and u^ ' are jointly Gaussian 
or that UjlUjU^'Uj’ satisfies the hypothesis of equation (12) since 


U-jUjU^Ui ' also enters in 


(u 




The Johns Hopkins University, 

Baltimore , Mi., June 9, 1952. 



38 


NACA TN 3116 


APPENDIX A 


Consider the fourth-order correlation u^UjU-^'u^' where is 

the velocity at the point, x and ujj.' is the velocity at x' = + JL* 

An outline of the procedure, based on invariant theory and essentially 
following Robertson (ref. 23), will be given for deriving the form for 
the isotropic tensor u^u^u^'^ ' . Consider the scalar correlation 

between the two velocity components in two arbitrary directions at x 
and two velocity components in two arbitrary directions at x ' . Let 

and b^ be the direction cosines of the two arbitrary directions at x 
and c^ and d^ be the corresponding quantities at x ' . Then the scalar 
correlation R is 


R(a,b,c,d;r) = u^UjUjj.'Uj * a^b^^c^d^ 


(Al) 


where r = (iplp) is the distance between two points x and x' . 

The correlation R(a,b,£,d;r) has the following special properties: 

(1) It is invariant under an arbitrary translation or rotation, as 

a rigid body, of the configuration defined by the points x and x ' and 
the unit vectors a, b, £, and d. 

(2) Its value is unchanged by the reflection of the above configu- 
ration in any point. These two are the conditions of homogeneity and 
isotropy. 

(3 ) It is homogeneous quadrilinear in the components of the four 
vectors a, b, £, and d. 

The form of R(a,b,£,d;r ) has to be determined under the above 
three conditions. According to the invariant theory of rotation groups 
in three dimensions (ref. 23), any invariant function of any number of 
vectors j[, a, b, . . . can be expressed in terms of the fundamental 
in variants of the following types: (l) The scalar product (_£ b) = 

of any two vectors including the scalar square and (2) the 

determinants 



£l a l *1 
^2 &2 1>2 
b a 5 b 
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of any three vectors. In terms of geometrical notions the invariants 
associated with a set of vectors are their lengths g.^, the angle 

between any two vectors, and the volume jj_ a hj of the par alle lepiped 
whose edges are any three given vectors. However, the volume of the 
parallelepiped changes sign on reflection, jj_ a b~] = - [a|_l|, hence 
these invariants do not appear in the expression for R(a,h,c,djr ). 

Since R(a,h,£,d;r) is a homogeneous quadrilinear in the components of 
the unit vectors a, b, £, and d, it must he a lin ear function of the 
form (| a)(i b)(j. c)(£ d), the six forms of type (l a)(_| h)(c d), and 

three of the type (a b)(£ d), with coefficients which are even functions 
of r. Thus, 

R (§^£>i> r ) = R i( r )(l £)(i 5.)(i c)(g d) + R 2 (r)(| a)(|_ b)(£ d) + 

£)(t d)(a h) + R^(r)(|_ d)(|_ b)(a £) + 

R 5 (r)(i a)(£ £)(b d) + Rg(r)(i b)(£ £ )(a d) + 

R y( r )(l a.)(£ d)(£ b) + Rg(r)(a b)(£ d) + 

R g(r)(a £)(b d) + Rio( r ^£ <*)(]> £.) 

~ R i( I * R 2^ 1 *^i^j^kZ a i^ ) j c k^'Z 
^(rJl^jh^aih^dj + R 4(r)l l i J 6 ik a i h (j c k d z + 

R 5 (r) ^i l k 5 jz a i b j c k d z + R 6( r ^^k 5 iz a i^ c k d i + 
RY(r)5 1 ^ 2 5j k a i hjC k d z + Rg(r jS^S^a^jC^ + 

E g( r )5ik5j Z a 1 hj c kd z + Rio( r ) 5 iz 5 kJ a i t ’j c k d Z (A 2 ) 

Equating equations (Al) and (A 2) and making use of the fact that this 
equality is true for arbitrary unit vectors a, b, c, and d. 
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'WV = + KaCrJl^ + V^WiJ + 

R Jj.( r ) & Z&j 5 lk + R 5( r )tilk 5 j2 + R 6^ r ^«J^k 8 i2 + 

Rj( r )iil z 8j z + R 8^ r ^ 5 ij S k2 + R 9^ r ^ 5 ik S ^l + R lo( r > 8 il8kJ 

(A3) 

This is ttie general form of the correlation involving four velocity com- 
ponents, either two velocity components at one point and two at another, 
or three velocity components at one point and one at another. Use is 

made of the fact that for the correlation u^u^u^'u^* the indexes i 

and j’ can he interchanged, k and 2 can he interchanged, and ij 
can he interchanged with k2. Using these symmetry conditions it is 
found that 


U i U j U k' U 2 I = ¥i 5 jM2 + K* 2 + R 3)(Mj 6 k2 + M2 & ij) + 

+ R 5 + ®6 + * 7 ) (Mj 8 ik + Mk 8 j2 + ^^k 5 i2 + 
5 i l 2 5 j2) + R 8 5 ij 6 k2 + |( R 9 + R lo)(&i2&kj + 


The five scalars characterizing u^uju^’u^ ’ can he expressed in terms 

of five special correlations. "When this is done equation (l6) is 
obtained. 


As mentioned earlier, equation ( 3 ) gives the form for the general 
fourth-order correlation, so that u^u^u^u^", the correlation involving 

three velocity components at one point and one at another point, has the 
form 


u i u j u k u 2 ' = S l( r ^i^j^2 + ^i^j 8 k2 + ^fr^k^ij + s 4^ r ^2^J 8 ik + 

S 5^ r ^i^k 6 j2 + S 6^ r ^j^ 8 i2 + S 7^ r ^i^2 5 jk + S 8^ r ^ 5 ij S k2 + 

s 9( r ) 6 U^j2 + S lo( r ^ 8 i2 8 k t 5 



Making use of the fact that Indexes i, j, and k can he interchanged. 


¥WV = S l(r)^i^Ml + + S 5^ + S 6( r J| (liSjB k2 + 

Mk 5 ,n + 6 j 6 k s n) + 3 [^ r ) + + (mz 6 ^ + 

Mj 5 ik + Mi 8 Jk) + 3 |® 8 ( r ) + s 9^ r ^ + B i0 ( r 3 ( 5 ij 5 kZ + 
6 ik B jZ + 5 iZ 5 kj) 


The four scalars defining u^UjU^u^ 1 can he expressed in terms of four 
special correlations. When this is done equation (19) is obtained. 
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TABLE I 

PRESSURE CORRELATIONS COMPUTED FROM FOURTH- 
AND SECOND- ORDER CORE® RATIONS 



P 2 / p2 ( U l 2 ) 2 

fef/'W 

(1) 


— //dp V 2 

s lv>k) 

(2) 

1/2 

Two extreme values 
computed from 
fourth-order 
correlation 

3-3 

.45 

85 

11 

0.28 

.20 

Computed from 
second-order 
correlation 

0.50 

12 

0.20 


dimension of values is approximately (in. ) 
^Dimension of values is approximately in. 












TABLE XI 

EXPERIMENTAL DATA 


Investigator and date 


In. 

in. 

2 A 2 

M 2 

25 , 2.3 

K* 2 R * 



\° y / 

IpVA 

m 

OJ 

Q. 

Simmons 1935 

57.0 

0.16 

0.08 

8.00 

0.08 

7.92 

0.35 

Collis 1948 

23.0 

.14 

.20 

1.00 

•15 

•85 

1.08 

Collie 1948 

38.0 

.14 

.20 

1.00 

.09 

•91 

1.05 

Referenc e ( 13 ) , 195 1 

29.0 

.17 

.11 

4.78 

.11 

4.67 

.46 

Do. 

32.0 

.12 

.12 

2.00 

.10 

1.90 

.73 

Do. 

36.O 

.41 

.36 

2.60 

.08 

2.52 

.63 

Do. — 

42.5 

.23 

.235 

1. 92 

.07 

1.9 

• 74 

Do. 

45.5 

.165 

.245 

.910 

.07 

.84 

1.09 

Do. 

49.0 

.55 

.22 

5.04 

.06 

4.98 

• 45 

Do. 

6l.O 

.23 

• 34 

.914 

.04 

.87 

1.06 

Do. 

67.O 

.34 

.22 

4.76 

.04 

4.72 

.46 

Do. 

74.0 

.19 

.27 

• 990 

.04 

.93 

1.02 
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